function [neighbors]= node_neighbors(node, matrix)

   neighbors = find(matrix(node,:) == 1); 
   temp=find(neighbors == node);
   neighbors(temp)=[];

end